First-order scaling near a second-order phase transition: 
Tricritical polymer collapse 
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The coil-globule transition of an isolated polymer has been well established to be a second-order 
phase transition described by a standard tricritical 0(0) field theory. We provide compelling evi- 
dence from Monte Carlo simulations in four dimensions, where mean-field theory should apply, that 
the approach to this (tri)critical point is dominated by the build-up of first-order-like singularities 
masking the second-order nature of the coil-globule transition: the distribution of the internal energy 
having two clear peaks that become more distinct and sharp as the tricritical point is approached. 
However, the distance between the peaks slowly decays to zero. The evidence shows that the posi- 
tion of this (pseudo) first-order transition is shifted by an amount from the tricritical point that is 
asymptotically much larger than the width of the transition region. We suggest an explanation for 
the apparently contradictory scaling predictions in the literature. 
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Mean-field theory is generally applicable to second- 
order phase transitions above their upper critical dimen- 
sion, and so is believed to provide an adequate descrip- 
tion of the approach to such critical points. One type 
of transition where mean-field theory should hold are tri- 
critical points |Q] for dimension d > 3. The region around 
a tricritical point in general dimension jjj is described by 
crossover scaling forms, where quantities depending on 
two relevant parameters can be essentially described by 
functions of a single scaling combination of those two pa- 
rameters. However, as noted by Lawrie and Sarbach (see 
page 106 of 0) there may be some breakdown of this 
crossover scaling for d > 3 due to the presence of one or 
more dangerous irrelevant variables. Here we provide ev- 
idence of a rather more dramatic breakdown of that sim- 
ple crossover scaling for the case of the coil-globule tran- 
sition of an isolated polymer, which is generally accepted 
to be described by tricritical theory We demon- 

strate that it is likely that the build up of the tricritical 
point is through the forming of singularities that have 
more in common with a (non-critical) first-order transi- 
tion! This however can be explained by a different kind 
of mean-field approach (not starting with an explicitly 
tricritical Landau functional) and, moreover, the region 
around the tricritical point needs to be described by more 
complex scaling forms. This second issue is in fact sepa- 
rate from the first-order nature of the scaling approach: 
we speculate that this behaviour is intimately related to 
the general description of systems where mean-field the- 
ory is used, so may have more general applicability. 

An isolated polymer in solution is usually argued to be 
in one of three states depending on the strength of the 



inter-monomer interactions which are mediated by the 
solvent molecules and can be controlled via the temper- 
ature T. At high temperatures and in so called "good 
solvents" a polymer chain is expected to be in a swollen 
phase (swollen coil) relative to a reference Gaussian state 
so that the average size R of the polymer scales with 
chain length faster than it would if it were behaving as 
a random walk. At low temperatures or in poor solvents 
the polymer is expected to be in a collapsed globular 
form with a macroscopic density inside the polymer: this 
implying an average size that scales slower [|| than a 
random walk, that is 
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so the globular state has radius-of-gyration exponent 
v g = 1/d. Between these two states there is expected to 
be a second-order phase transition (sharp in the infinite 
chain length limit). The standard description of the col- 
lapse transition is that of a tricritical point related to the 
n — > limit of the 4> 4 -(p 6 O(n) field theory . One then 
might expect that above the upper critical dimension 
(d u — 3) some type of self-consistent mean-field theory 
based upon a suitable tricritical Landau-Ginzberg Hamil- 
tonian would give a full description of the transition, 
and hence conclude that in all dimensions d > 3 there is 
a collapse transition from a swollen state to the globular 
state with classical tricritical behaviour. In the corre- 
spondence P] between a Landau-type functional for a 
magnetic system with a magnetisation M of n compo- 
nents (in the limit n — * 0) and the polymer problem, the 
coefficient of the M 4 in the Landau functional maps to 
the second virial coefficient B of the polymer solution and 



1 



so to the temperature T of the polymer, while the tem- 
perature of the magnet r is related to the polymer length 
N through N ~ (r — r c ) _1 where t c (B) is the location of 
the phase transition line in the so-called symmetry plane 
jlj of the tricritical system. A finite polymer length im- 
plies the corresponding magnet is at a high temperature 
r > t c . Hence, the finite-size scaling around the collapse 
transition should be described by tricritical crossover. 

The application of the mean-field theory of a tricritical 
point to polymer collapse predicts that at the transition 
point the polymer actually behaves as if it were a ran- 
dom walk (vg — 1/2), and this point has been known 
as the 0-point. Thermodynamically (N = oo), one ex- 
pects a weak transition with a jump in the specific heat 
a = (note that the thermodynamic polymer exponent 
a is related to the shift exponent in tricritical theory [jj] , 
itself not to be confused with the polymer theory finite- N 
scaling ip p shift exponent). For finite polymer length N 
there is no sharp transition for an isolated polymer (un- 
less one examines a macroscopic number of such poly- 
mers) and so this mean-field transition is rounded and 
shifted. In three dimensions the application of various 
self-consistent mean-field like approaches predicts that 
the second-order transition is rounded and shifted on the 
same scale of N~ x / 2 , that is, the crossover exponent <fi 
is 1/2, though strictly the power laws involved are mod- 
ified via renormalisation group arguments by confluent 
logarithms. In four and higher dimensions no confluent 
logarithms should be present and one may expect pure 
mean-field behaviour with a crossover exponent of 1/2 || 
(cj>t is the analogous tricritical exponent p]]). 

On the other hand, for d > 3 Sokal [0 has pointed out 
that the alternative method of analysing collapse which 
has been shown to be equivalent to the field theoretic ap- 
proach, namely the continuum Edwards model, has dif- 
ficulties: in fact, if one analyses the Edwards model one 
finds the crossover exponent is given by <pE — 2 — d/2, 
which for d = 4 gives (f>E = 0! In passing we note here 
that the same analysis predicts the shift of the 0-point, 
defined say via the universal ratio of the radius of gy- 
ration to the end-to-end distance equalling its Gaussian 
value, should scale as aM^ 2 " 1 ) so ip E = (d/2-1) ^ <t> E - 
This difference between the shift and the crossover ex- 
ponent implies that strict crossover scaling has broken 
down. Of course, the theoretical fact that the swollen 
phase should also be Gaussian for d > 4 does raise the 
suspicion that the analysis of the Edwards model for 
polymer collapse may be subtle for d > 3. 

To resolve this question of the crossover scaling one 
may first be tempted to try some simple scaling argu- 
ments as follows. Around the 6*-point (T — » Tg and 
N — > oo) a crossover form 0] is generally predicted to 
be: 

R N ~ N 1/2 T{{T e - T)N*). (2) 
The mathematical understanding of these forms was re- 



cently re-examined 0. To match the behaviour of ([!]) 
consider the behaviour of @ for T < Tg and fixed, so we 
need to find the behaviour of T(z) for z — > oo. Assum- 
ing a power law T(z) ~ z , and the mean-field density 
exponent (3 = 1, one finds that b = 1/d and more im- 
portantly that <fr = (d/2 — 1). So we notice that the 
crossover exponent found from this simple scaling argu- 
ment is the same as predicted from the Edwards model 
for the shift exponent. Hence we might conclude that 
for d > 3 the crossover scaling form might still be appli- 
cable but with shift and crossover exponents that obey 
<\> = %jj p = (d/2 — 1). While this implies a subtlety in 
the Edwards model analysis and also one in the map- 
ping of the tricritical theory to the polymer problem, it 
is perhaps not completely surprising. 
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FIG. 1. Internal energy density distributions for lengths 
N = 2048 and 16384, each at their respective transition tem- 
peratures. The more highly peaked distribution is associated 
with length 16384. 



To consider such issues we have simulated interacting 
self-avoiding walks, the canonical lattice model of poly- 
mer collapse, on the four-dimensional hyper-cubic lattice 
using the PERM algorithm [|| over a wide range of tem- 
peratures with surprising results. Because this algorithm 
is based upon kinetic growth it works well around the col- 
lapse region as well as throughout the swollen phase. As 
a consequence we are able to obtain reasonable data up 
to length N — 16384 (a more complete discussion of our 
simulations is available in |^]). Now, firstly, our results 
suggest that there is indeed a collapse transition in four 
dimensions at a finite temperature. However, the charac- 
ter of that transition is particularly intriguing! In d = 4 
our scaling argument above predicts <p — ip p — 1 while 
the mapping to the tricritical theory predicts <p = 1/2: 
we find neither! In fact we find a rounded transition 
with a divergent specific heat, and near the transition the 
distribution of internal energy is clearly bimodal. This 
'double-peaked' distribution becomes more pronounced 
as the chain length is increased (see figure 1). As we vary 
temperature through the transition region both peaks 
are essentially stationary and one grows as the other 
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decreases in size: classic first-order behaviour. On the 
other hand we also were able to find a candidate 6>-point 
(a critical state) where R ~ N 1 / 2 well above the transi- 
tion region. We then considered the shift of the rounded 
first-order like transition to the 6>-point: the best scaling 
produced a shift exponent of about 1/3 (see figure 2). 
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FIG. 2. Scaling of the shift of the transition: we show the 
scaling combination iV 1 ' 3 (Wc,jv — <^e) versus N~ 2 ^ 3 , where 
Uc,n and u)e are the Boltzmann weights associated with the 
monomer-monomer interaction at T Ci at and Tg. 

By studying the position when the universal ratio of 
the mean square distance of a monomer to the end-point 
to the mean square end-to-end distance takes on its Gaus- 
sian value we find that the 0-point is shifted much less 
and may scale as 1/N. Hence there seem to be two shift 
exponents! While these results seem at variance with 
standard tricritical ideas there is a mean-field type the- 
ory that describes the first-order transition region well. 
This framework was explained some time ago in a paper 
by Khokhlov jl0|], who applied the mean- field approach 
of Lifshitz, Grosberg and Khokhlov O [lj| to arbitrary 
dimensions. Here we argue that the conclusions of these 
works may be valid for d > 4. The theory is based on a 
phenomenological free energy in which the competition 
between a bulk free energy of a dense globule and its 
surface tension drive the transition. The consequences 
of this surface free energy were largely ignored in the 
polymer literature until recently, when its effect on the 
scaling form of the finite-size partition function was pro- 
posed and confirmed |i^-|l7||. 

Lifshitz theory JllJ is based on several phenomenologi- 
cal mean-field assumptions. Firstly, there exists a 0-state. 
Secondly, for lower temperatures there exists a globular 
state where the polymer behaves as a liquid drop. The 
results of the theory are based on a phenomenological 
free energy of that globular state relative to the free en- 
ergy of the pure Gaussian state of the 0-point at Tg. The 
starting point of this analysis is a bulk free energy with a 
quadratic dependence on the distance to the 0-point and 
so an exponent a = 0, implying a second-order phase 
transition occurs in the thermodynamic limit. From the 



theory one finds a rounded transition for finite N occur- 
ring at temperature T c ^ < Tg with shift exponent 

^„ = l/(d-l) (3) 

so that the collapse occurs at finite length at a temper- 
ature that scales towards the 0-point quite slowly and is 
below the #-point. This then concurs with our finding 
that ip p « 1/3 in d — 4. The width of this transition AT 
at finite N can be found and has crossover exponent 

0=(d-2)/(d-l). (4) 

All the exponents are derived from the assumptions of 
mean field thermodynamic behaviour and of eq. (jl]). 
Hence all the exponents quoted here are related to each 
other (only one independent exponent) . Our simulations 
show just such a narrow crossover region with a crossover 
exponent that is certainly larger than the shift exponent 
(the data is compatible with exponent 2/3 in d = 4). 

By considering the difference of the density p s of the 
swollen state to the globular state p g relative to the den- 
sity of the swollen state at T Ci jv, Khokhlov [[To) concluded 
from its divergence that 'the coil-globule transition is 
first-order', though we now interpret this to mean that 
the finite-size corrections to the thermodynamic second- 
order transition are first-order like. We point out that 
the terminology of Khokhlov was presumably that ex- 
plained in Section I.C.2 of jl3| but may be misleading to 
the modern reader. However, both p g {T Ct iq) and p s (T c ^) 
tend to zero as N — > oo and it is simply that p g {T c ^) 
tends to zero asymptotically slower than p s (T c ^) that 
makes the relative difference diverge. The analysis can 
be used to deduce the scaling of Rn at T c at with an ex- 
ponent v c = l/(d— 1). Note that vg > v c > v g , so that 
this scaling is in-between the scaling fixed at the 0-point 
and at any temperature fixed in the collapsed phase. Fol- 
lowing the work of Lifshitz, Grosberg and Khokhlov (ll| 
one can also calculate the change in the internal energy 
over the crossover width of the transition AT as the la- 
tent heat AU from the free energy expression. The la- 
tent heat decays as N increases with exponent 1/(1 — d). 
The corresponding height of the peak in the specific heat 
diverges with exponent (d — 3)/(d— 1). From our sim- 
ulational data we were unfortunately unable to extract 
reasonable estimates for v c or the exponents of the latent 
and specific heats. 

To interpret our results we can take the understand- 
ing of this mean-field theory further. Let us consider 
the distribution of internal energy, which we measured in 
our simulations, as a function of temperature and length. 
For any temperature above the 0-point and well below 
T Cj tv one expects the distribution of internal energy to 
look like a single peaked distribution centred close to the 
thermodynamic limit value: a Gaussian distribution is 
expected around the peak with variance OiN^ 1 ! 2 ). In 
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fact this picture should be valid for all temperatures out- 
side a range of order 0(iV~( d ~ 2 )/( d_1 )) centred on T c jv. 
When we enter this region we will expect to see a dou- 
ble peaked distribution as in a first-order transition re- 
gion. For any temperature in this transition region there 
should be two peaks in the internal energy distribution 
separated at the order of (^(iV -1 /^ -1 )) (the value of the 
gap being the latent heat). Each of these peaks should be 
of Gaussian type with individual variances again of the 
order 0{N~ 1 ^ 2 ). Hence as N increases the peaks will 
become more and more distinct and relatively sharper 
but the peak positions will be getting closer together. 
Hence we refer to this scenario as a psewdo-first-order 
transition or, more correctly, as first-order-like finite-size 
corrections to a second-order phase transition. If there 
were a real first-order transition then the distance be- 
tween the peaks should converge to a non-zero constant. 
A comprehensive interpretation of our computer simu- 
lations in d = 4 jj^l is most consistent with just such a 
scenario and leads us to conjecture that this theoretical 
picture is indeed correct for the coil-globule transition for 
d > 4. 

Let us return to the question of crossover scaling forms 
and our finding of two shift exponents from our computer 
simulations. While we cannot ascertain either with great 
accuracy let us assume that we have a region around the 
#-point that is approximated well by a form like (||) with 
crossover exponent 4>e = (<i/2 — 1) . Now, despite the fact 
that this does not describe the collapse transition region, 
we notice that substituting t ~ iV 1 /^ -1 ), and using the 
asymptotics derived from our matching argument above, 
leads to R ~ JVVld-i^ which is precisely the correct 
scaling for the real transition region! Hence we conjecture 
a phenomenological product scaling form (for T < Tg) 

R N ~ N 1 ' 2 F((T e - T)N^)g{{T^ N - T)N+) (5) 

with T C . N ~ Tg - aN-^ and where Q(y) ~ 1 for 
y — > ±oo. This form will then correctly describe both 
the region around the 0-point and the rounded transi- 
tion around T Ci jy and will match with the behaviour of 
the collapsed phase for fixed T <Tg. Such a form is not 
dependent on the finding of pseudo-first-order behaviour 
and may be useful for analysing data whenever two shift 
exponents are found. Two shift exponents may appear 
in systems that are described by mean-field theories. 

Finally let us re-examine the crossover exponent de- 
rived from the Edwards model. This is essentially only 
valid for T > Tg. For d > 5 the swollen phase is Gaussian 
itself and the polymer temperature is an irrelevant vari- 
able (<j)E < 0), and in d = 4 it is marginal, so one might 
suspect that the scaling theory for T > Tg is very differ- 
ent from that for T < Tg. Hence while the scaling form 
(^) from the Edwards model with (j> — 4>e correctly pre- 
dicts that R ~ iV 1 / 2 /(T) for d > 5 it does not describe 
the collapse. 



In conclusion, Monte Carlo simulations of lattice poly- 
mers in four dimensions show that for finite length the 
rounded coil-globule transition appears first order but we 
argue that the tricritical predictions may well reappear in 
the infinite length limit as our results fit best the predic- 
tions of Lifshitz-Grosberg-Khokhlov (LGK) theory ap- 
plied to high dimensions. We suggest that the crossover 
scaling forms are more complicated than at low dimen- 
sions and suggest a generalisation that may be heuris- 
tically useful. The Edwards model is valid above the 
0-point while the coil-globule transition is described by 
LGK theory. 
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